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ABSTRACT 

Aims. The detection of old neutron stars and stellar mass black holes in isolation is one of the cornerstones of compact object 
astrophysics. Microlensing surveys may help on this purpose since the lensing mechanism is independent of the emission properties 
of the lens. Indeed, several black hole candidates deriving through microlensing observations have been reported in the literature. The 
identification of counterparts, especially in the X-rays, would be a strong argument in favor of the compact nature of these lenses. 
Methods. We perform a cross-correlation between the catalogs of microlensing events by the OGLE, MACHO and MOA teams, 
and those of X-rays sources from XMM-Newton and Chandra satellites. Based on our previous work, we select only microlensing 
events with duration longer than one hundred days, which should contain a large fraction of lenses as compact objects. Our matching 
criterion takes into account the positional coincidence in the sky. 

Results. We find a single match between a microlensing event, OGLE 2004-BLG-81 (t E ~ 103 days), and the X-ray source 2XMM 
J180540.5-273427. The angular separation is ~ 0.5 arc-seconds, i.e. well inside the 90% error box of the X-ray source. The hardness 
ratios reported in the 2XMM catalog would suggest a hard spectrum with a peak between 2 and 4.5 keV or a softer but highly absorbed 
source. Moreover the microlensing event is not fully constrained, and possible association of the source star with a flaring cataclysmic 
variable or a RS Canum Venaticorum-like star have been suggested as well. 

Conclusions. The very small angular separation (within uncertainties) is a strong indicator that 2XMM J 180540.5-273427 is the 
X-ray counterpart of the OGLE event. However, the uncertainties on both the nature of lensed system and the lens itself challenge 
the interpretation of 2XMM J180540.5-273427 as the first confirmed isolated black hole identified so far. In any case, it probes the 
diagnostic capacity of microlensing surveys and open the path for further identifications of black hole or neutron star candidates. 
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1. INTRODUCTION 

The Galactic population of neutron stars and and black holes 
(NSs and BHs h ereafter) is expected to c ount as much as 10 8 - 
10 9 objects (e.g. lKeane & Kra mer. 2008). This population is al- 
most undetected, with the exception of few thousands of sources, 
mainly isolated young NSs emitting as radio/X-ray/gamma-ray 
pulsars and NSs and BHs in X-ray binaries. 

Observing the large majority of old, likely isolated objects 
is difficult since the proposed recycling mechan ism, i.e. resid- 
ual accretion from the interst ellar medium (e.g. lOstriker et all 
[197(1 iTreves & ColpiL 1 19911) , is hampered by many factors , 
like the large velocities ( e.g. ITreves et all l2000t IPopov et"aTJ 



120001: ISartore et all 1201 , t he pre sence of the mag netic field 
in NSs (e.g. lToropinaetall 120031: iPerna et all 120031) and, for 



BHs, the small radiation efficiencies of the accretion proces s 
(e.g. lAgol & Kamionkowskii l2002t iBeskin & Karpovl 120051) . 
Nevertheless, isolated old remnants (NSs and BHs) can certainly 
act as gravitational lenses. 

Gravitational microlensing (ML) events occur when light 
from a star is de flected and amplified by a foreground mass 
(|Paczvnskiill986l) . Indeed, several thousan ds of ML events have 
been observed so far (see iMoniezL 1201 Ol for a review). Since 
most of the events are due to stars with mass less than one so- 
lar mass, NSs and BHs contribute essentially to the lo ng duration 
event s , the longest ones bein g associated to B Hs (e.g.lMao et all 
I2002t IBennett et all 120021) . According to iBennett et all the 
number of long duration events is in excess of that inferred 



from the oretical models (e.g. | Gouldl l2000HWood & Maoll200l . 
Recently Sartore & Treves! (120 10t) showed that this excess can be 
explained assuming that BHs receive large natal kic ks at birth, of 
the same magnitude of those observed in NSs (e.g. lHobbs et all 
l2005tlGualandris et all 120051) . Under this assumption, ~ 40 per- 
cent of the events longer than 100 days are due to BHs. 

Thus, microlensing surveys may have already detected tens 
of events related to NSs and BHs. However, the nature of can- 
didate BH lenses has to be assessed with independent methods. 
The mass estimated from ML events is in general affected by 
large uncertainties, given the degeneracy of the time-scale with 
the input parameters like the mass of the lens and the geome- 
try of the source-lens-observer system, i.e. the relative positions 
and velocities. On the other hand, the detection of X-rays com- 
ing from the position of the ML event would be a strong hint 
for the presence of a compact object. A search for the X-ray 
counterpart of one of the most promising BH candida t es de- 
tected t hrough ML has been perform ed by iNucita et al.l |2006, 
see also Agol & Kamionkowski 2002). However, a 100 ks point- 
ing with XMM-Newton gave no result down to a limiting flux 
~ 10~ 15 erg s in the 0.2 - 10 keV energy range. 

In this work we present a systematic cross-correlation analy- 
sis of long duration ML events with the X-ray catalogs of XMM- 
Newton and Chandra satellites, which appeared recently. The 
paper is organized as follows: Sect|2] shortly describes the ex- 
pected spectra of NSs and BHs accreting at low rates. In Sect|3] 
we report on the list of ML events and catalogs of X-ray sources 
used to find possible counterparts. A description of the cross- 
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correlation method is also given. In Sect 4 we present our re- 
sults and we briefly discuss on a source resulting form the cross- 
correlation procedure. Finally in Section [5] we summarize the 
main points of the paper and discuss the possibilities of finding 
more remnant candidates with ML. 

2. EXPECTED X-RAY SPECTRA OF ACCRETING 
REMNANTS 

It is commonly accepted that isolated NSs and BHs accret- 
ing from the ISM should radiate mostly in the X-ray band. 
For instance, the spectrum emitted by a NS accreting at low 
rate could be crudely approximated to a Planckian, with a 
temperature of f ew hundreds elettronvo l t (lOstriker et all Il970t 
iTreves & ColpiL 1 199 it iBlaes & MadauL 1 19931) More detailed 
calculations dZampieri et al.L 119951; IZane et all 120001) showed 
that the spectrum emerging from an accreting NS, either magne- 
tized or non-magnetized, could have substantial deviations from 
a blackbody and in particular the resulting spectrum would be 
harder than the blackbody emission at the same accretion rate. 
However, the overall behavior would still resemble a Planckian- 
like spectrum with peak emission in the soft X-rays. 

The situation for BHs is quite different since, at variance 
with NSs, the efficiency of conversion of accreted matter into ra- 
diation is expected to be lower. Models of ISM accretion onto 
isolated BHs have been develo p ed by lAgol & K amion kowskil 
(120021) : I Beskin & Karpovl (l2005): Map elli et alJ(l2006h . The un- 
derlying physics is assumed to be similar to that of BH can- 
didates in X-ray binaries but the exact details are poorly con- 
strained. Nevertheless, the emerging spectra should have to have 
both thermal-like and non thermal components, as observed in 
X-ray binaries. 

The X-ray luminosities of accreting NSs/BHs would depend 
heavily on the relative velocity of the comp act object (NS /BH) 
with respect to the surrounding medium (e.g. lBondilll952l) . 
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where 77 is the efficiency of the X-ray emission, n is the num- 
ber density of the ISM, v and M are the relative velocity with 
respect to the surrounding medium and the mass of the compact 
objects, respectively. Old, but still fast moving NSs are expected 
to po pulate the Galactic disk with velocities of ~ 200 km s" 1 
(e.g. lSartore et all uOlO). The resulting luminosity is of the order 
of ~ 10 27 - 10 28 erg s _I . For BHs, the effect of their higher mass 
is counterbalanced by the lower efficiency and thus the nominal 
X-ray luminosities are approximatel y of the same order of ac- 
creting NSs, see however the paper of Mape lli et all Also, when 
comparing with observations one has to take into account the 
absorption effect, which are severe in the soft X-rays. 

3. METHOD 

3.1. Selection of microlensing events 



We ba se our work on the public data o f the OGLE ( Udalski et al 



1992), MACHO (lAlcock et all Il99l and MOA (Murakietal 
1999) collaborations available on the World Wide Web 



Following the results of Sarto re & TrevesL we select only events 
with time-scale longer than 100 days. The OGLE (Optical 
Gravitational Lensing Experiment) datcQ were collected with the 



Early Warning System (EWS, see lUdalskiL 120031) from 1998 
to 2009 and correspond to the OGLE-II and OGLE-III phases 
of the survey. The number of events detected in the first year 
(1998) was 41, of which none had a duration longer than 100 
days while e.g. in 2008 the number of events detected was 654 
with 38 events longer than 100 days. The total number of events 
from OGLE is 4117, of which 177 percent fall in our range of 
interest. The MACHO (MAssive Compact Halo Objects) sur- 
vey data3 were collected from 1993 to 1 999 and are compre- 
hensive of 528 bulge events reported by Tho mas et all d2005l) 
plus th e Red Clump Giants events reported by Popowsk i et all 
(2005), totaling 567 events. Of these, 69 are longer than 100 
days. The MOA (Microlensing Observations in Astrophysics) 
surve)Q started in 2000, and in 2006 in entered in its second 
phase. While the number of events detected in 2000 was barely 
8, the current detection rate is similar to that of the OGLE-III 
survey, i.e. ~ 500 - 600 yr _I . The number of events detected by 
MOA up to 2010 is 2622, of which 268 show time-scales longer 
than 100 days. 

3.2. X-ray sources catalogs 

We searched for X-ray counterparts of our selected sample of 
ML events in the version 1.2 of the second XMM-Newton 
Seren dipitous Source Catalog (2XMM hereafter, IWatson et all 
2009) and in t he version 1 . 1 of the Chandra Source Catalog 
(CSC hereafter, Evans et al., 2010). The most recent version of 
the XMM catalogFjhas been released in April 2010. It contains 
data of 191870 unique sources with a median flux (in the 0.2 - 
12 keV band) ~ 2.5 x 10~ 14 erg s _1 cm -2 , with 20 percent of the 
sources having fluxes below 10~ 14 ergs~' cirT 2 . The CSC cata- 
log released in August 2010, covers approximately 320 square 
degrees of the sky at the 10~ 13 erg s cirT 2 flux limit (0.5 - 7.0 
keV band). The sky coverage dr ops to ~ 6 square d egrees for 
a flux limit of 10~ 15 erg s -1 cirT 2 ( lEvans et aLl TJoiO). The total 
number of unique sources contained in the catalog is 106586. 

The OGLE-III survey monitored a region of the sky com- 
prised between I = -2° and I = -7° and b = ±10°, i.e. ~ 100 
square degrees. The area covered by the MACH O and MOA-II 
surve ys is roughly half of that of OGLE-III. (e.g. lThomas et all 
2005). The XMM-Newton and Chandra surveys made respec- 
tively 17 and 5 pointings towards the same are£0 (see Fig. [TJ 
which, considering the 30' x 30' field of view of the EPIC cam- 
eras on board XMM-Newton and the 16' x 16' field of view of 
the ACIS cameras on board Chandra, correspond to ~ 14 and 
~ 1 square degrees, respectively. 

3.3. Cross-correlation analysis 

Given the large uncertainties affecting the expected observa- 
tional appearance of old isolated NSs and BHs, we do not apply 
any particular spectral or variability criteria on possible X-ray 
counterparts. Nevertheless, we get rid of spurious sources in the 
2XMM catalog by selecting only those with S UM.FLAG = 0, 
i.e. those for which none of the warning flags of the EPIC instru- 
ments were set. These flags indicate the probability of a detected 
source to be spurious. For the CSC catalog there is no analogue 



2 http://wwwmacho. mcmaster. ca/ 



http://www.phys. canterbury, ac. nz/moaj 



http://ogle. astrouw. edu.pl/. 



4 http://xmmssc-www. star. le. ac. uk/Catalogue/xcat public J2XMM. html 

5 http://cxc.harvard.edu/csc/ 

6 multiple observations of the same region are considered as a single 
pointing 
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flag, but the estimated number of spurious sources is expected 
to be less than one in every field with 100 kiloseconds of inte- 
gration. Thus, our matching criterion is based only on the po- 
sitional coincidence of the long time-scale ML events with the 
X-ray sources of the 2XMM and CSC catalogs. 

Our cross-correlation software computes the projected dis- 
tance between the selected ML events and the entries of both the 
2XMM and CSC catalogs. A positive match is found when a ML 
event lies with the 3<x error circle of an X-ray source. Positional 
errors of single sources are taken from the respective catalogs. 
The positional error of 2XMM sources (POS ERR column) al- 
ready accounts for systematic errors. The total lcr uncertainty is 
calculated as t he squared root of th e sum of statistical and sys- 
tematic errors (Wat son et all 12009) 



POS ERR = ^iRADEC-ERR 2 + S YSERR 2 . 



(2) 



For CSC sources, to take into account the systematic error 
(0.16 arc-seconds) we use the equation suggested by the CSC 
team 



err.ellip.r 0ttot = 2.4477467 x (3) 
x ^0.1669041 x {e r r.ellip.n Xcat ) 2 + 0.0256 . 



Typically, the statistical error for on-axis CSC sources is ~ 
0.2 arc-seconds, while for off-axis sources at 14 arc-minutes the 
statistical error is ~ 3.5 arc-seconds. We assume that positional 
errors of ML events are of cr ML ~ 1 .5 arc-seconds for all the 
events. Thus, the resulting radius of the error circle is assumed 
to be the root mean square of the X-ray source and ML event 
positional uncertainties 



ML 



+ crt 



(4) 



4. RESULTS 

The cross-correlation analysis returned a single positive match in 
the 2XMM catalog. The associated lensing event was observed 
in 2004 by both the OGLE and MOA surveys and is identified as 
OGLE 2004-BLG-81 and MOA 2004-BLG-3, respectively. The 
duration of the event reported by the OGLE team is ~ 103.63 
days. However, the li ght curve is poorly fitted b y standard lens- 
ing models (Fig. |2j. Wvrzvkows ki et all J2006) found that the 
baseline of the source star, i.e. the magnitude outside the ML 
event (I ~ 17), has a suspected periodicity of ~ 4 days, thus 
pointing to an eclipsing binary. In particular, the shape of the 
folded light curve points to a contact binary system. 

To add more confusion, the MOA team reports a baseline 
/ ~ 8 (sic!), a duration of ~ 6.73 days and amplification very 
close to unity, A ~ 1.002. However, a visual inspection of the 
stellar field does not confirm the presence of such a bright star, 
whose image would have been affected by diffraction. Thus we 
rely exclusively on the OGLE data to characterize the event. 

The X-ray source associated to the ML event, 2XMM 
J 180540.5-273427 (J 1805 hereafter), has been serendipitously 
observed during a pointing of MACHO-96 -BLG-5, another BH 
candidate detected through microlensing dBennett et all [2002; 



iNucita et all 12006)) . The X-ray properties of the source have 
been retrieved with the XCat-DB web interface dMotch et all 
2009). The total number of counts is 312.744 + 0.001 
(0.2 - 12 keV band) corresponding to a flux of (3.39 ± 
0.78) x 10~ 14 erg s _1 crrT 2 , which implies a luminosity, ne- 
glecting the photoelectric absorption of the ISM, of ~ 3 x 
10 30 {d/1 kpc) 2 erg s _1 . If J 1085 is the responsible for the magni- 
fication of a bulge star (d ~ 8 kpc), then it should be placed at an 
intermediate distance and its X-ray luminosity should be lower 
than ~ 10 32 ergs~'. 

The positional uncertainty of the source is ~ 2 arc-seconds 
and it lies at ~ 0.5 arc-seconds from the position of the ML 
event. We report the fluxes on the different EPIC bands and the 
relative hardness ratios as given in the XCat database in TableQ] 

5. DISCUSSION AND CONCLUSIONS 

The small angular separation (~ 0.5 arc-seconds) between the 
position of the ML event OGLE 2004-BLG-81 and J1805 is 
well below the lcr positional uncertainty of the X-ray source and 
makes the association highly likely. Thus, if J 1805 is actually a 
BH, it would prove that ML surveys can detect isolated compact 
objects. However, there are a number of uncertainties that need 
to be addressed in order to not discard the claim. 

First, the nature of the event reported by the OGLE and MOA 
surveys is unclear. As already pointed out, the shape of the light 
curve would rule out a ML event, even accounting for secondary 
effects like parallax or blending. The fact that the source star 
is possibly a contact binary would indicate a cataclysmic vari- 
able (CV), i.e. an accreting white dwarf. This would imply that 
the event was in reality an outburst episode rather than a gen- 
uine gravitational lens. An accreting white dwarf would easily 
explain the detected X-ray radioation as coming from the mat- 
ter accreted by the degenerate star. However, the duration of the 
putative outburs t and it s lightcurve are unusual for this kind of 
sources (see e.g. lButlertll99lHKuulkers et alll2003l) . thus chal- 
lenging the CV hypothesis. 

In alternative, it has been sugg ested that the len sed star is a 
chromospherically active variable (Ber nhardt 120091) . possibly of 
the RS Canum Venaticorum (RS CVn) type. This class of vari- 
ables is known to show periodic variations which are thought to 
be related to the active regions on the surface of the star or ellip- 
ticity of the star itself. These effects can mimic the lightcurve of 
an eclipsing binary and explain the observed modulation of the 
baseline. Furthermore, RS CVn stars are known X-ray emitters, 
with luminosities of ~ 10 31 ergs~', and they also show flaring 
activity at both optical and X-ray wavelengths. Yet, also in this 
case the amplitude, shape and duration of the optical event are 
unusual for RS CVn stars (S.N. Shore, private communication). 

Aside from the possible classification of the optical event, 
the nature of the X-ray source J 1805 is poorly constrained as 
well. The small number of photons collected does not allow a 
good characterization of the spectrum. The hardness ratios are 
positive at low energies and change sign between 2 and 4.5 keV. 
This fact would suggest a hard or highly absorbed spectrum. A 
hard spectrum would rule out a NS as accreting object since its 
temperature is expected to be below ~ 1 keV. Thus, if J 1805 is 
an isolated compact object, it should be a BH. 

In conclusion, in this paper we performed a cross-correlation 
analysis of long time-scale microlensing events with the soft X- 
ray sources of the recently released XMM-Newton and Chandra 
catalogs. We report a positive match between the ML event 
OGLE 2004-BLG-81 and the X-ray source 2XMM J180540.5- 
273427. The nature of both the ML event and the the X-ray 
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Table 1. Observed properties of the X-ray source 2XMM J 180540.5-273427. The data are retrieved from the XCat database. 
Hardness ratios are calculated as described in 1 Watson et all with the lower energy band corresponding to the position in the table. 



Energy Band 


0.2-0.5 


0.5-1.0 


1.0-2.0 


2.0-4.5 


4.5-12 


[keV] 












Flux 


0.002 ± 0.001 


0.034 ± 0.027 


0.298 ± 0.050 


1.20 ±0.146 


1.73 ±0.762 


[xl0" 14 ergs~' cnT 2 ] 












Hardness ratio 


0.982 ± 0.268 


0.738 ±0.141 


0.132 ±0.100 


-0.356 ±0.157 





source are still unclear. Follow-up optical and X-ray observa- 
tions of the field of 2XMM J 180540.5-273427 at high angular 
resolution could help to characterize the source star and the X- 
ray source as well. 

We stress that, even if the association between OGLE 2004- 
BLG-81 and 2XMM J180540.5-273427 would not be con- 
firmed, the argument that a large fraction of long duration events 
are related to BHs is still valid. As shown in Section [2] and in 
particular in Equation [TJ the expected flux from an isolated BH 
or NS depends on many unknown quantities: the velocity of the 
collapsed object, the density of the ISM, the photoelectric ab- 
sorption etc. Therefore, it is not straightforward to use the ab- 
sence, or paucity, of correlation between X-ray sources and ML 
events as a constraint on so many parameters. 

In the next years we expect a substantial enlargement of 
the catalogs of ML events and X-ray sources. In particular the 
eROSITA mission, which should be launched shortly, will make 
a survey of the entire sky, which in the soft X-ray band (0.5 
- 2 keV) will be 30 times more sensitive than ROSAT. At the 
same time, deep systematic exposures with Chandra and XMM- 
Newton, of the most promising ML events, chosen mainly on the 
basis of the duration, will also help to set stronger constraints on 
the flux of accreting isolated compact objects. 
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Fig. 2. Light curve of the event OGLE 2004-BLG-8 1 . Solid line 
represents the best fit to the light curve assuming a standard mi- 
crolensing model. Source: OGLE website. 



Fig. 1. Centroids of the pointings, in Galactic coordinates (I, b), 
made by the XMM-Newton (top) and Chandra (bottom) satel- 
lites in the region of the sky monitored by the OGLE-III survey 
(see text). 
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